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Abstract.  We  evaluate  the  effects  of  solid  particle  rotation  in  high-altitude  solid  rocket  exhaust  plume  flows,  through  the 
development  and  application  of  methods  for  the  simulation  of  two  phase  flows  involving  small  rotating  particles  and  a 
nonequilibrium  gas.  Green’s  functions  are  derived  for  the  force,  moment,  and  heat  transfer  rate  to  a  rotating  solid  sphere 
within  a  locally  free-molecular  gas,  and  integration  over  a  Maxwellian  gas  velocity  distribution  is  used  to  determine  the 
influence  of  particle  rotation  on  the  heat  transfer  rate  at  the  equilibrium  limit.  The  use  of  these  Green’s  functions  for  the 
determination  of  particle  phase  properties  through  the  Direct  Simulation  Monte  Carlo  method  is  discussed,  and  a 
procedure  is  outlined  for  the  stochastic  modeling  of  interphase  collisions.  As  a  test  case,  we  consider  the  nearfield  plume 
flow  for  a  Star-27  solid  rocket  motor  exhausting  into  a  vacuum,  and  vary  particle  angular  velocities  at  the  nozzle  exit 
plane  in  order  to  evaluate  the  influence  of  particle  rotation  on  various  flow  properties.  Simulation  results  show  that 
rotation  may  lead  to  slightly  higher  particle  temperatures  near  the  central  axis,  but  for  the  case  considered  the  effects  of 
particle  rotation  are  generally  found  to  be  negligible. 


INTRODUCTION 

In  the  modeling  and  analysis  of  exhaust  plumes  from  solid-propellant  rockets,  several  complex  and  poorly 
understood  physical  phenomena  are  often  overlooked.  In  particular,  the  AI2O3  particles  which  result  from  the 
combustion  of  aluminized  propellants  are  generally  considered  using  simplified  two  phase  flow  models,  due  to  a 
lack  of  experimental  study  or  numerical  methodologies  for  some  of  the  more  complicated  physical  processes 
involved.  One  such  example  is  particle  rotation;  little  study  has  been  devoted  to  the  effects  of  particle  rotation  in  the 
flow  regimes  common  to  rocket  exhaust  plumes,  and  rotation  has  generally  been  neglected  in  plume  modeling  and 
analysis.  In  order  to  understand  and  quantify  the  influence  of  particle  rotation  on  various  flow  properties,  the 
mechanisms  by  which  particles  in  the  plume  develop  angular  momentum  must  first  be  considered. 

In  a  typical  solid  rocket  motor,  aluminum  is  used  within  the  propellant  grain  to  reduce  combustion  instabilities 
and  increase  specific  impulse.  Liquid  alumina  (ALO3)  droplets  are  formed  inside  the  combustion  chamber,  with  a 
bimodal  size  distribution  concentrated  around  diameters  of  roughly  1  pm  and  100  pm  [1],  The  droplets  are  then 
accelerated  through  the  rocket  nozzle,  where  they  may  encounter  several  different  phenomena  that  tend  to  contribute 
angular  momentum.  First,  smaller  droplets  are  accelerated  by  the  expanding  gas  far  more  rapidly  than  the  larger 
droplets,  so  that  droplets  of  different  sizes  develop  widely  varying  velocities.  Alumina  droplets  commonly  account 
for  20%  to  40%  of  the  total  mass  flow  through  the  nozzle  entrance,  so  the  high  number  density  and  considerable 
relative  velocities  between  droplets  of  different  sizes  allow  for  a  large  number  of  collisions.  Most  of  these  collisions 
result  in  coalescence  of  the  droplets  into  larger  agglomerates,  which  rotate  according  to  the  dynamics  of  off-center 
coalescing  collisions  [2].  At  the  same  time,  larger  droplets  tend  to  break  apart  within  the  nozzle  due  to  collisions  or 
surface  stresses  imposed  by  the  accelerating  gas,  and  these  breakup  processes  may  impose  significant  angular 
momentum  on  the  resulting  droplet  products.  Within  the  subsonic  converging  nozzle  region,  larger  droplets  may 
diverge  considerably  from  gas  streamlines  and  accumulate  in  the  boundary  layers  along  the  nozzle  wall,  where  large 
transverse  gradients  in  the  bulk  gas  velocity  can  also  induce  rotation  [3],  As  the  droplets  cool  and  solidify  into  solid 
particles,  inelastic  collisions  between  particles  may  contribute  to  rotation  as  well.  Turbulent  and  acoustic 
fluctuations  in  the  gas  can  intensify  many  of  these  rotation  mechanisms. 
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In  exhaust  plume  flows  at  high  altitudes,  the  gas  generally  exhibits  a  high  degree  of  nonequilibrium  behavior 
outside  of  a  small  inviscid  core  region  just  beyond  the  nozzle  exit  [4].  The  gas  velocity  and  energy  distributions  can 
deviate  significantly  from  the  equilibrium  limits,  and  continuum  assumptions  underlying  Computational  Fluid 
Dynamics  (CFD)  techniques  may  break  down  due  to  the  large  mean  free  paths  in  relation  to  the  characteristic  length 
scales  for  macroscopic  flow  gradients.  In  order  to  accurately  model  such  flows,  the  nonequilibrium  nature  of  the  gas 
must  therefore  be  considered  within  a  non-CFD  framework.  As  an  additional  problem  in  the  consideration  of  two 
phase  flows,  typical  CFD-based  methods  for  two  phase  flow  simulation  use  rarefaction  corrections  to  empirical 
particle  drag  coefficient  and  Nusselt  number  formulations  which  often  fail  at  the  high  Knudsen  number  and  high 
Mach  number  conditions  of  interest. 

As  an  alternative  to  empirical  correlations  for  interphase  momentum  and  energy  transfer,  analytical  models  for 
solid  bodies  in  high  Knudsen  number  flows  may  provide  significant  improvements  in  both  overall  accuracy  and  the 
physical  effects  which  may  be  considered.  Particle  rotation  is  one  such  effect,  and  several  authors  have  studied 
various  properties  associated  with  the  rotation  of  objects  within  a  free-molecular  equilibrium  gas.  Epstein  [5] 
derived  the  moment  imparted  on  a  rotating  sphere,  subject  to  several  constraints  including  a  lack  of  translational 
motion.  Wang  [6]  determined  the  drag  and  lift  coefficients  for  a  rotating  sphere,  and  showed  that  the  direction  of  the 
lift  force  is  opposite  that  of  the  Magnus  force  for  continuum  flow.  Ivanov  and  Yanshin  [7]  derived  force  and 
moment  equations  for  various  bodies  of  revolution.  While  all  of  these  authors  have  developed  closed  form  solutions 
which  are  applicable  to  the  general  flow  regimes  of  interest,  the  formulas  they  provide  are  limited  by  assumptions  of 
Maxwellian  gas  velocity  distributions,  and  none  have  considered  the  effect  of  rotation  on  interphase  heat  transfer.  In 
addition,  these  formulas  tend  to  be  very  complicated  and  poorly  suited  for  evaluation  within  a  larger  numerical 
scheme  for  the  simulation  of  high-altitude  plume  flows. 

In  a  recent  paper  [8]  a  methodology  is  proposed  for  the  simulation  of  such  plume  flows,  where  the  gas  is 
modeled  using  the  Direct  Simulation  Monte  Carlo  (DSMC)  method  and  the  rates  of  momentum  and  heat  transfer  to 
the  particles  are  based  on  equations  of  Gallis  et  al.  [9]  While  the  particles  are  assumed  to  be  chemically  inert, 
consideration  is  made  for  the  effects  of  interphase  momentum  and  energy  transfer  on  both  the  particles  and  the  gas. 
In  the  present  effort,  this  method  is  extended  to  include  effects  of  particle  rotation,  so  that  the  influence  of  rotation 
on  various  flow  properties  can  be  evaluated.  First,  we  derive  the  Green’s  functions  for  the  force,  moment  and  rate  of 
heat  transfer  on  a  rotating  spherical  particle  in  a  locally  free-molecular  gas.  For  the  special  case  of  a  stationary 
rotating  particle  in  a  Maxwellian  gas,  the  heat  transfer  Green’s  function  is  integrated  analytically  to  provide  a  closed 
form  solution.  Next,  we  discuss  the  use  of  these  Green’s  functions  to  modify  the  properties  of  representative  solid 
particles  within  a  two-phase  DSMC  simulation,  and  present  two-way  coupling  approaches  through  which  solid 
particles  may  influence  the  surrounding  gas.  When  used  within  a  larger  two-phase  flow  model,  these  procedures 
allow  for  momentum  and  energy  conservation  in  a  time-averaged  sense,  and  enable  the  simulation  of  two-phase 
rarefied  flows  involving  very  large  solid  particle  mass  fractions. 

The  procedures  mentioned  above  are  implemented  for  two  phase  flow  calculations  in  the  DSMC  code 
MONACO  [10].  As  a  test  case,  we  consider  the  near-field  plume  flow  for  a  Star-27  solid  rocket  motor  [11] 
exhausting  into  a  vacuum.  Estimates  of  particle  angular  velocities  at  the  nozzle  exit  are  based  on  the  assumption  that 
coalescing  droplet  collisions  within  the  nozzle  are  the  dominant  mechanism  for  inducing  particle  rotation,  and 
calculations  are  made  using  the  limiting  cases  of  maximum  possible  and  zero  initial  angular  velocities.  Comparison 
of  results  shows  that,  at  least  for  the  case  considered,  effects  associated  with  particle  rotation  are  generally  within 
the  limits  of  statistical  scatter  in  the  calculated  flow  properties. 


ANALYSIS  OF  ROTATION  EFFECTS 

Following  the  Green’s  function  approach  of  Gallis  et  al.  [9],  we  consider  a  spherical  solid  particle  within  a 
locally  free-molecular  simple  gas.  Incident  gas  molecules  which  collide  with  the  particle  are  either  specularly 
reflected  off  the  surface,  or  are  diffusely  reflected  with  full  thermal  accommodation  of  translational  and  rotational 
energy  modes  to  the  particle  temperature.  We  assume  that  the  particle  is  small  enough  to  allow  for  a  low-Biot 
number  approximation  of  spatially  uniform  temperature,  yet  large  enough  that  the  change  in  particle  velocity  during 
an  interphase  collision  is  much  smaller  than  that  of  the  gas  molecule.  These  assumptions  are  all  generally  valid  for 
the  high  altitude  plume  flows  of  interest,  and  are  used  to  derive  the  following  Green’s  function  for  the  force  on  a 
nonrotating  particle  due  to  gas  molecules  of  relative  velocity  class  ur: 

F„>r)  =7rR^n/(ur )fmcr  +^27imkBTp  (!) 
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Here  Rp  is  the  particle  radius,  ng  is  the  gas  number  density,  /(ur)  is  the  gas  velocity  distribution  function  in  a  particle- 
centered  coordinate  frame,  m  is  the  mass  per  gas  molecule,  cr  is  the  relative  speed,  x  is  the  particle  thermal 
accommodation  coefficient,  Tp  is  the  particle  temperature  and  kB  is  Boltzmann’s  constant. 

In  considering  the  Green’s  function  Fp(ur)  for  the  force  imparted  by  the  gas  on  a  rotating  particle,  we  first 
decompose  Fp  into  a  term  which  accounts  for  momentum  transfer  to  the  particle  from  incident  gas  molecules,  a  term 
for  momentum  exchange  during  specular  reflection,  and  a  corresponding  term  for  diffuse  reflection.  It  can  be  shown 
that  the  second  term  is  zero  and  that  only  the  third  term  may  be  influenced  by  particle  rotation.  Next,  using  the 
velocity  superposition  principle,  the  momentum  transfer  due  to  diffuse  reflection  is  separated  into  a  component 
independent  of  rotation  and  a  component  which  depends  on  the  particle  angular  velocity  cop  but  is  independent  of  Tp. 
This  latter  component  is  the  only  contributor  in  a  time-averaged  sense  to  the  force  Green’s  function  Frot  due  to 
particle  rotation,  and  allows  the  Green’s  function  for  the  total  force  on  the  particle  to  be  computed  as  Fp  =  Fnr  +  Frot. 
Following  the  above  logic,  Frot  is  equal  to  the  product  of  the  collision  frequency  for  diffusely  reflecting  collisions 
and  the  rotational  component  -m<Ut>  of  the  average  momentum  transferred  to  the  particle  during  diffuse  reflection. 
Here  ut  is  the  tangential  velocity  of  the  particle  surface  at  the  collision  point  in  an  inertial  particle-centered  reference 
frame.  In  order  to  evaluate  <ut>  we  express  ut  as  a  function  of  the  angle  9  between  -ur  and  the  outward  surface 
normal  unit  vector  nc  at  the  collision  point,  then  multiply  by  the  distribution  function  for  9  and  integrate  over  all 
particle  surface  elements  where  a  collision  may  occur.  We  find: 

2n 

Fro,(ur)=yRjmn/(ur)xcc1,x^  (2) 


By  adding  (1)  and  (2)  and  integrating  over  all  gas  velocity  space,  we  then  determine  the  total  force  on  the  particle. 

In  the  absence  of  any  gas  vorticity  effects  or  forces  not  associated  with  interphase  collisions,  the  rotating  particle 
will  experience  a  damping  moment  about  the  particle  center  which  will  tend  over  time  to  reduce  the  magnitude  of 
CQp.  Depending  on  the  angle  between  cop  and  the  gas  bulk  velocity  relative  to  the  particle,  the  moment  may  also  have 
a  component  orthogonal  to  cop  which  tends  to  change  the  axis  about  which  the  particle  rotates.  This  moment  can  be 
determined  using  the  same  Green’s  function  approach  as  above.  Again,  incident  and  specularly  reflected  gas 
molecules  do  not  contribute  to  the  moment,  and  the  average  moment  contribution  during  diffuse  reflection  depends 
only  on  the  tangential  velocity  of  the  particle  surface.  The  Green’s  function  Mp(ur)  for  this  moment  can  then  be 
expressed  as  the  negative  product  of  the  collision  frequency  for  diffuse  reflection  and  the  average  angular 
momentum  mRp<ncxut>  about  the  particle  center  imparted  on  diffusely  reflected  gas  molecules.  The  averaging 
operation  is  performed  through  integration  over  the  particle  surface,  to  find  the  following  result: 


^ur)=^mn/(ur)xcr 


1^-3® 


(3) 


V  ~r  J 

Through  energy  conservation  arguments,  the  Green’s  function  for  the  heat  transfer  rate  to  a  rotating  particle  can 
be  expressed  as 


QP(u,)=  ^Rp  n/Klxc, 


1  2 

2mC:-^n+l 


(4) 


where  erot  is  the  average  rotational  energy  of  incident  gas  molecules,  Arot  is  the  number  of  rotational  degrees  of 
freedom  for  the  gas,  and  ekin  =  2kBTp+  !/2m<ut-ut>  is  the  average  kinetic  energy  of  diffusely  reflected  gas  molecules 
in  an  inertial  particle-fixed  coordinate  frame.  By  integrating  over  the  particle  surface,  we  derive  the  general  formula 

/  o\1 

^mc;-(  2  +  kBi;  +erot  +  R; 

^  '  \  1  v  yj 

where  ®p=|(Op|.  Note  from  (5)  that  a  greater  value  of  (op  corresponds  to  an  increase  in  heat  transfer  from  the  gas  to  the 
particle.  This  property  can  be  justified  physically  by  considering  a  single  particle  surface  element  in  a  coordinate 
frame  which  is  fixed  to  the  local  tangential  surface  velocity.  In  this  coordinate  frame  the  average  kinetic  energy  of 
diffusely  reflected  gas  molecules  will  be  independent  of  a>p,  while  the  average  kinetic  energy  of  incident  gas 
molecules  tends  to  increase  as  a>p2.  Thus,  the  only  effect  of  particle  rotation  here  is  to  increase  the  interphase  energy 
transfer  associated  with  the  kinetic  energy  of  incident  gas  molecules  which  are  involved  in  diffusely  reflecting 
collisions. 

In  the  special  case  where  the  particle  which  moves  at  the  bulk  velocity  of  an  equilibrium  simple  gas  with 
temperature  Tg,  (5)  may  be  integrated  analytically  over  all  gas  velocity  space  to  provide  a  closed  form  solution  for 
the  net  heat  transfer  rate: 


QpK)=^Rjn/(RK 


Approved  for  public  release;  distribution  unlimited. 


qj,rct=^ngT^27ikBTg/m  (4+Arot)kB(Tg-Tp)+|ni£i|R; 


If  Tp  >  Tg  as  is  typical  for  high-altitude  rocket  exhaust  plume  flows,  then  the  two  bracketed  terms  in  (6)  are  of 
opposite  sign.  Assuming  the  first  term  is  greater  in  magnitude  than  the  second,  we  find  that  particle  rotation  tends  to 
reduce  the  magnitude  of  the  heat  transfer  rate  and  allows  for  a  more  gradual  reduction  in  Tp  over  time. 


DSMC  IMPLEMENTATION 


The  Green’s  functions  given  above  can  be  used  to  allow  rotating  solid  particles  to  be  modeled  within  a  DSMC 
simulation.  For  simplicity  we  neglect  radiative  heat  transfer  and  assume  the  particle  phase  to  be  dilute  and 
chemically  inert,  so  that  temporal  variation  in  the  properties  of  individual  particles  may  occur  only  through  the 
mechanisms  described  above  or  through  particle-wall  interactions.  Representative  solid  particles  are  tracked  along 
with  DSMC  gas  molecules  through  a  computational  grid,  and  during  each  time  step  the  total  force,  moment,  and 
heat  transfer  rate  for  each  particle  is  determined  by  adding  contributions  from  all  gas  molecules  which  are  assigned 
to  the  same  cell.  The  force,  moment  and  heat  transfer  contributions  to  a  particle  from  a  single  DSMC  gas  molecule 
are  computed  by  evaluating  (1),  (2),  (3)  and  (5)  with  the  quantity  ng/(ur)  replaced  by  the  ratio  of  the  local  gas 
molecule  weighting  factor  to  the  cell  volume.  Once  all  Ng  molecules  in  the  cell  have  been  considered,  the  particle 
velocity  up  and  temperature  Tp  are  altered  as  described  in  Ref.  [8],  and  cop  is  updated  by  the  product  of  the  net 
moment,  the  local  time  step  and  the  inverse  of  the  particle  moment  of  inertia. 

Due  to  the  large  particle  phase  mass  fractions  typical  of  the  high-altitude  plume  flows  of  interest,  accurate 
simulation  requires  consideration  of  momentum  and  energy  coupling  from  the  particles  to  the  gas.  This  is 
accomplished  here  by  probabilistically  modeling  individual  interphase  collisions  and  modifying  the  velocity  and 
rotational  energy  of  all  gas  molecules  involved.  A  collision  selection  scheme  based  on  the  No  Time  Counter  method 
of  Bird  [12]  is  used  to  determine  which,  if  any,  DSMC  gas  molecules  will  collide  with  each  solid  particle  during  a 
given  time  step.  If  a  particular  gas  molecule  is  found  to  collide  with  the  particle,  then  the  collision  will  involve  either 
diffuse  reflection,  with  probability  x,  or  specular  reflection,  with  probability  1-x.  In  the  latter  case,  the  post-collision 
relative  velocity  u,  is  sampled  from  an  isotropic  distribution  such  that  the  relative  speed  cr  of  the  molecule  is 
unchanged  during  the  collision,  and  the  molecule  velocity  is  reassigned  to  equal  up+ur  .  In  the  case  of  diffuse 
reflection,  a  more  complicated  procedure  is  required  to  find  ur ,  and  the  gas  molecule  rotational  energy  must  be 
altered  as  well. 

If  a  gas  molecule  is  diffusely  reflected  off  a  nonrotating  solid  particle,  then  the  magnitude  and  direction  of  ur  are 
statistically  independent,  so  each  can  be  determined  separately.  Following  the  procedure  outlined  in  Ref.  [8],  the 
post-collision  relative  speed  cr  is  sampled  from  the  distribution  function 

/(c;)dc;=2(^c^exp(-^c;2)dc;  (7) 

using  the  acceptance-rejection  method  [12],  where  (3p  =^m/(2kBT ) •  The  angle  5  between  the  vectors  -ur  and  u,  is  also 

determined  using  the  acceptance-rejection  method,  by  sampling  from  a  numerical  approximation 

/(5)  =  0.02  0  4  256  -  0.25  1  555+  1.10454-  1.90  3  53  +  0.493  852+  1.24  85  (8) 

of  the  distribution  function  for  5e[0,  ji].  An  azimuthal  angle  a0  required  to  define  the  direction  of  ur  in  three 
dimensional  space  is  sampled  from  a  uniform  distribution  over  [0,  2ji],  and  the  components  of  ur  are  calculated 
using  formulas  derived  by  Bird  [12]  for  binary  elastic  collisions: 

v,-4-v,«B8-gg5g^iB5S)l.  ^J^+"SgSggg3>]  (9) 

'  c,l  ■  '  J  c,l  tf+ufr  J  <4  J 

The  absolute  velocity  of  the  gas  molecule  is  then  set  to  the  final  value  up+ur  . 

If  a  diffusely  reflecting  collision  involves  a  rotating  solid  particle,  then  Eq.  8  is  not  valid  and  an  alternate 
procedure  must  be  used.  In  this  case  u,  is  found  by  separately  determining  nonrotational  unr  and  rotational  urot 
components,  then  adding  the  two  together.  To  calculate  unr ,  we  first  compute  the  angle  9  between  -u,  and  the 
outward  normal  unit  vector  nc  at  the  collision  point  on  the  particle  surface  by  setting  sin20  equal  to  a  random  number 
on  [0,  1].  A  corresponding  azimuthal  angle  ai  is  sampled  from  a  uniform  distribution  over  [0,  2ji],  and  components 
of  nc  are  then  calculated  as 

+  n  =i[  n  =l[-wrC069+  ^v^q-u^sin^l  (i0) 

crL  J  ?  Cr|_  J  cr|_  (v[+wr)/! 
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Next,  we  find  the  angle  C,  between  the  vectors  nc  and  unr  by  sampling  sin2C  from  a  uniform  distribution  over  [0,  1], 
An  azimuthal  angle  a2  is  sampled  from  a  uniform  distribution  over  [0,  2ji],  and  Eq.  7  is  used  to  find  the  magnitude 
of  unr  .  The  components  of  unr  are  then  determined  through  the  following  equations: 


j  ’  (l-x/  j  "  d-t,r  j 

Finally,  the  rotational  component  of  u,  is  computed  as  urot  =  Rp  copxnc ,  and  the  gas  molecule  velocity  is  reassigned 
to  the  post-collision  value  up+  unr  +  urot  . 


PLUME  SIMULATION 

To  evaluate  the  importance  of  particle  rotation  effects,  all  procedures  described  above  are  implemented  in  the 
DSMC  code  MONACO  [10],  and  an  axisymmetric  simulation  is  performed  for  the  nearfield  plume  flow  from  a 
Star-27  solid  rocket  motor  expelling  into  a  vacuum.  The  nozzle  has  an  exit  diameter  of  78  cm  and  a  divergence  half¬ 
angle  of  17.2°.  All  gas  and  particle  properties  along  the  nozzle  exit  plane  are  taken  from  Anfimov  et  al.  [11]  Based 
on  values  provided  by  these  authors,  spherical  A1203  particles  of  diameter  Dp  =  0.15,  0.2,  0.3,  0.5,  1,  2  and  3  pm 
together  constitute  a  total  mass  flow  fraction  of  about  30%,  with  particles  of  each  size  distributed  uniformly  over  the 
nozzle  exit.  The  gas  here  is  a  mixture  of  H2,  N2  and  CO,  with  mole  fractions  of  0.38,  0.31,  and  0.31  respectively. 
Effects  of  radiative  heat  transfer,  heterogeneous  reactions,  gas  chemistry  and  particle  phase  change  are  neglected. 
While  each  of  these  phenomena  may  significantly  alter  flow  properties  of  interest,  we  intend  only  to  assess  the 
importance  of  particle  rotation  effects  on  a  representative  plume  flow,  not  to  calculate  flow  properties  for  a 
particular  case  with  the  greatest  possible  accuracy. 

Following  Vasenin  et  al.  [2],  we  assume  that  coalescing  collisions  between  liquid  droplets  within  the  nozzle  are 
the  dominant  contributors  to  particle  rotation.  The  maximum  allowable  particle  angular  velocity  magnitude  is  then 
determined  by  a  criterion  for  the  centrifugal  breakup  of  liquid  droplets.  A  normalized  angular  velocity 


is  set  to  equal  some  critical  value  Qcr;t  for  all  particles  at  the  nozzle  exit,  where  mp  is  the  particle  mass  and  a  is  the 
liquid  surface  tension  of  the  particle  material.  While  the  magnitude  of  cnp  will  be  constant  for  all  particles  of  equal 
diameter  at  the  nozzle  exit,  the  initial  direction  of  cop  for  each  particle  is  randomly  selected  on  a  plane  normal  to  the 
axis  of  symmetry.  Salita  [1]  uses  momentum  and  energy  conservation  arguments  to  show  that  the  collision-induced 
rotation  rate  of  spherical  liquid  agglomerates  should  be  limited  by  Qcrit=3.31.  Following  these  arguments,  we 
perform  plume  simulations  with  initial  particle  angular  velocities  corresponding  to  Q=3.31  and  Q-()  in  order  to 
evaluate  potential  upper  bounds  for  the  influence  of  particle  rotation  on  various  flow  properties.  For  both  cases, 
particle-to-gas  coupling  is  considered  through  application  of  the  method  described  above  involving  (10)  and  (11),  to 
avoid  any  differences  in  results  caused  by  the  numerical  approximation  in  (8). 

The  simulations  use  a  rectangular  grid  domain,  starting  at  the  nozzle  exit  plane  and  extending  10  m  downstream 
in  the  axial  direction  and  4  m  in  the  radial  direction.  Calculations  are  performed  on  16  1.4  GFlz  AMD  Opteron 
processors,  with  each  case  requiring  about  3000  CPU  hours.  As  an  example  of  simulation  results  for  the  Q=3.31 
case,  Figure  1  shows  average  angular  velocity  magnitudes  for  particles  of  each  size  as  a  function  of  distance  from 
the  central  axis,  along  a  radial  plane  10  m  downstream  of  the  nozzle  exit.  Figure  2  shows  the  variation  in  average 
particle  temperatures  along  this  same  plane,  for  both  the  Q=3.31  and  Q=0  cases.  From  a  comparison  of  TP  values  for 
any  radial  location  and  particle  size,  particle  rotation  is  shown  in  Figure  2  to  have  little  measurable  impact  on 
particle  temperatures  in  the  plume.  Flowever,  for  all  but  the  largest  and  smallest  particles  considered,  rotation  does 
account  for  a  roughly  5  to  10  K  increase  in  average  particle  temperatures  within  about  0.5  m  of  the  central  axis.  This 
trend  is  likely  due  to  the  fact  that  the  gas  density,  and  therefore  the  component  of  the  interphase  energy  transfer  rate 
associated  with  particle  rotation,  tends  to  decrease  with  distance  from  the  axis.  The  trend  is  less  pronounced  for  the 
smallest  particles  because  these  rapidly  lose  the  bulk  of  their  rotational  energy  on  entering  the  grid  domain,  and  for 
the  largest  particles  because  of  their  relatively  low  initial  angular  velocities. 

With  this  one  minor  exception,  we  find  a  general  lack  of  dependence  on  particle  rotation  for  all  flow  properties 
of  interest.  Distributions  of  gas  translational  temperature,  particle  and  gas  number  densities,  maximum  divergence 
angles  for  each  particle  size,  and  various  other  properties  show  differences  between  the  two  cases  which  are 
comparable  to  the  expected  statistical  scatter.  We  can  therefore  conclude  that,  at  least  for  the  plume  simulation 
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considered  here,  particle  rotation  may  be  neglected  with  no  significant  loss  in  overall  accuracy.  The  numerical 
methods  described  above  do  maintain  some  utility  in  making  this  conclusion  possible,  and  in  allowing  for  the 
simulation  of  other  two-phase  rarefied  flows,  such  as  MEMS  flows  involving  frequent  particle-wall  collisions, 
where  rotation  effects  may  still  be  important. 
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FIGURE  1.  Variation  in  particle  angular  velocity  magnitudes  along  a  radial  plane  10  m  downstream  from  the  nozzle  exit. 


:^*****************±*±±i^i***^*^m§JS 

’  '5OoOOOOOOOOOOOOOOOOOOOOOO  aAaA 

aaaa&aaa^aaA 

.  iAAAAAAAAAAAA AAAAAA&AA 


1.E+03-  □□□ 


□□□□□□□□□□□nan0000 


□  □□ 


□  □ 


□  □□ 


>0000<>0<>000^<><>0' 


o^' 


❖❖<> 


oo° 


ooO 


❖O' 


❖O 


1800 

1600 

1400 
„  1200 
^  1000 
800 

600 

400 

0  12  3  4 

y(m) 


****: 


*** 


** 


:  ^**yq*++*i  i  *  $  *  q-m4****** 


****** 


**** 


lWK.WWW*’M’" 

;  AAAAA  A  AAAAAAAAAAAAAAAAAAAAAAAAA^AA 


Dp  Gwm),  Q 

*  0.15,3.31 
O  0.15,0 

■  0.2,3.31 
□  0.2,0 

*  0.3,3.31 
A  0.3,0 

*  0.5,3.31 
o  0.5.0 

□  1,3.31 
x  1,0 
0  2,3.31 
+  2.0 
-3,3.31 
x  3,0 


FIGURE  2.  Average  particle  temperature  as  a  function  of  distance  from  the  central  axis. 
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